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Optical microrheology using rotating laser-trapped particles 

Alexis I. BishopQ Timo A. Nieminen, Norman R. Heckenberg, and Halina Rubinsztcin-Dunlop 
Centre for Biophotonics and Laser Science, Department of Physics, 
The University of Queensland, Brisbane QLD 4072, Australia 

We demonstrate an optical system that can apply and accurately measure the torque exerted by 
the trapping beam on a rotating birefringent probe particle. This allows the viscosity and surface 
effects within liquid media to be measured quantitatively on a micron-size scale using a trapped 
rotating spherical probe particle. We use the system to measure the viscosity inside a prototype 
cellular structure. 
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PACS numbers: 87.80.Cc,83.85.Cg,83.85.Ei 

It is well known that light can transport and transfer 
angular momentum as well as linear momentum Q]. In 
recent years, this has been widely applied to rotate mi- 
croparticles in optical tweezers. Such rotational micro- 
manipulation has been achieved using absorbing parti- 
cles j2,0,^irefringent particles 0,01, specially fabricated 
particles G,'^, and non-spheric al p articles using linearly 
polarized l ight 0, 0, 0, 0, 0, or non-axisymmetric 
beams ■ Notably, particles rotated in rotationally 

symmetric beams spin at a constant speed, indi- 

cating that the optical torque driving the particle is bal- 
anced by viscous drag due to the surrounding medium. It 
has been noted that the optical torque can be measured 
optically |l6l |. allowing direct independent measurement 
of the drag torque. 

For a pure Gaussian beam that is circularly polarised, 
each photon also has an angular momentum of ±h about 
the beam axis, equivalent to an angular momentum flux 
of ztP/w for the beam, where P is the power of an in- 
cident beam, and lo is the angular frequency of the op- 
tical field, which can be transferred to an object that 
absorbs the photon, or changes the polarization state 
The reaction torque on a transparent object that changes 
the degree of circular polarization of the incident light is 
given by 



TR = AaP/uj 



(1) 



where Aa is the change in circular polarisation. Thus, 
by monitoring the change in circular polarization of light 
passing through an object, the reaction torque on the 
object is found. We exploit this simple principle as the 
basis for an optically driven viscometer. The viscosity 
of a liquid can be determined by measuring the torque 
required to rotate a sphere immersed in the liquid at 
a constant angular velocity — a concept we have imple- 
mented at the micrometer scale using a system based on 



optical tweezers. We have rotated birefringent spheres 
of synthetically-grown vaterite that were trapped three- 
dimensionally in a circularly polarized optical tweezers 
trap and measured the frequency of rotation, and the 
polarization change of light passing through the particle, 
to determine the viscosity of the surrounding fluid. 

At least four distinct techniques have been used previ- 
ously for microrheology: tracking the diffusion of tracer 
particles 0,0], tracking the rotational diffusion of disk- 
shaped tracer particles |l9l | , measurement of correlations 
in the Brownian motion of two optically trapped probe 
particles [i^, and the measurement of viscous drag act- 
ing on an optically trapped probe particle in linear mo- 
tion |23|. The method we use here — direct optical mea- 
surement of the rotation of a spherical probe particle in 
a stationary position — allows highly localized measure- 
ments to be made, since the probe particle does not move 
in the surrounding medium, and, compared with the ro- 
tational diffusion measurements by Cheng and Mason us- 
ing microdisks 19], the use of spherical probe particles 
greatly simplifies the theoretical analysis of the fluid flow. 
Also, since the probe particle is rotationally driven by the 
optical trap, the rotation rate can be readily controlled. 

The expression for the drag torque, td, on a sphere 
rotatin g in a fluid for low Reynolds number flows is well 
known |22l | , and is given by 



(2) 



where /i is the viscosity of the fluid, a is the radius of the 
sphere, and Q is the angular frequency of the rotation. 
Equating the expressions for torques tr and tjj gives: 



fi = AaP/{8TTa^nuj). 



(3) 



This equation suggests that the viscosity can be de- 
termined using knowledge of only a few elementary 
parameters — the power and the polarization change of 
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the incident beam passing through the particle, and the 
size and rotation rate of the particle. The power, po- 
larization state, and rotation rate can all be measured 
optically. 

A significant obstacle to implementing this approach 
has been obtaining spherical, transparent, birefringcnt 
particles of suitable size for trapping (approximately 1~ 
10 /im in diameter). We have developed a novel tech- 
nique to grow nearly perfectly spherical crystals of the 
calcium carbonate mineral vaterite, which has similar 
birefringence properties to calcite, within the required 
size regime. We produce the vaterite crystals by adding 
6 drops of 0.1 M K2CO3 to a solution of 1.5 ml of O.IM 
CaCla plus 4 drops of O.IM MgS04. The solution is 
strongly agitated by pipetting. The solution initially ap- 
pears milky, becoming clearer as the crystals form within 
a few minutes of mixing. The typical mean size is found 
to be 3 niR with a spread of approximately ±1.5 /im. Fig- 
ure^shows both an optical image and an electron micro- 
graph of typical crystals used for viscosity measurements. 
Vaterite is a positive uniaxial birefringcnt material, with 
TT-o — 1.55 and rie — 1.65. 

(a) (b) 




FIG. 1: (a) Optical microscope image of a typical vaterite 
crystal used for viscosity measurements, (b) Scanning elec- 
tron microscope image of a vaterite crystal. 

A typical optical tweezers arrangement was used as the 
basis for the experiment, which delivered approximately 
300 mW of light at 1064 nm into a focal spot with a di- 
ameter of around 0.85 /zm using a lOOx oil-immersion 
objective of high numerical aperture (NA = 1.3). A A/4 
plate located immediately before the objective converts 
the linearly polarized light from the laser into a circu- 
larly polarized beam which is able to rotate birefringcnt 
objects 0. Light passing through the probe particle is 
collected by an oil coupled condenser lens (NA — 1.4) 
and then sent to a series of polarization measuring de- 
vices. The degree of circular polarization remaining in 
the beam is determined by passing it through a A/4 plate 
with the fast axis orientated at 45° to the axes of a po- 
larizing beamsplitter cube. Photodiodes monitoring the 
orthogonal outputs of the beamsplitter provide a mea- 
sure of the degree of circular polarization remaining in 
the beam, which allows the change in angular momen- 



tum, and hence the torque apphed to the particle to be 
known. The photodiodes are calibrated such that they 
provide a measure of the power at the trap focus. The 
sum of the two signals gives the total trapping power, 
and the difference between the signals gives the degree 
of circular polarization of the beam measured in units of 
power at the trap focus. The torque applied to the parti- 
cle is determined by measuring the difference between the 
initial degree of circular polarization of the beam in the 
absence of the particle and the final polarization, in ac- 
cordance with equation . The frequency of rotation of 
the particle can be accurately determined by monitoring 
the transmitted light through a linear polarizer, as sug- 
gested by Nieminen et al. [l6j. A small amount of light 
is diverted from the main beam and is passed through a 
polarizing beamsplitter cube, and the forward direction 
is monitored using a photodiode. The signal is modu- 
lated at twice the rotation rate of the particle, and the 
depth of the modulation is proportional to the amount 
of linearly polarized light in the beam. If the particle is 
thick enough to act as a A/2 plate, then circularly po- 
larized light passing through the particle is reversed in 
handedness and the frequency of the rotation cannot be 
measured as there is no linear component to the trans- 
mitted light. However, as the particles are spherical and 
are of similar diameter to the beam, the polarization at 
different radial distances from the rotation axis varies, 
and there will almost always be some linear component 
after transmission. 

The calculated viscosity depends on the cube of the 
radius of the particle, and thus it is important to deter- 
mine the size of the particles accurately. The particles 
used for viscosity measurements are in the range of 1.5 to 
3.5 /im in diameter and consequently to obtain the vis- 
cosity within 10% requires the diameter to be measured 
with an accuracy on the order of 90 nm, and it was found 
that obtaining an accurate measurement of the diame- 
ter by direct visualization was problematic. The method 
of measurement preferred by the authors is to put two 
spheres of nearly identical size in contact, and measure 
the distance between centres, which can be found ac- 
curately due to the spherical symmetry. The diameter 
is inferred using elementary geometry, and the error in 
the measurement by this method is estimated as being 
±40 nm. The degree of asphericity for typical particles 
was estimated to be less than 3%. 

The viscosity of water was measured using the system 
described. A range of vaterite particles were suspended 
in distilled water and trapped three dimensionally in a 
sealed slide/coverslip cell, with a depth of 50 /im. Parti- 
cles that are trapped begin to rotate immediately. How- 
ever, the rotation may be stopped by aligning the A/4 
plate located before the objective so that the polariza- 
tion is made linear. The size of the particles ranged from 
1.5 to 3.5 fim in diameter, and rotation rates up to 400 Hz 
(24000 rpm) were observed for powers up to 350 mW at 
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the trap focus. The Reynolds number of the fluid flow 
around the rotating spheres is quite low, on the order of 
1 X 10^'^ for rotations up to 1000 Hz, and hence the flow 
is well within the creeping flow regime required by 
The signals from the circular polarization and linear po- 
larization monitoring photodiodes were recorded using a 
16 bit analogue to digital converter sampling at 10 kHz 
for a period of 5 seconds. Typical signals recorded dur- 
ing an experiment are reproduced in figureEl The strong 
uniform modulation of the linear polarization signal (fig- 
ure |2Ia)) shows that the particle rotates very uniformly, 
taking a large fraction of the angular momentum carried 
by the beam. Using a particle of diameter 2.41 /zm, the 
viscosity was found to be /i = (9.6±0.9) x 10~'*Pas, which 
is in excellent agreement with the established value for 
the viscosity of water at 23°C of 9.325 x lO'^'Pas A 
series of measurements was made with a range of different 
sized particles at a range of rotation rates and trapping 
powers. The variation in viscosity over a range of powers 
below 100 mW was found to be on the order of 0.7% us- 
ing a single particle. The variation in inferred viscosity 
using 11 particles, spanning a range of sizes, was approx- 
imately 4.7%. The most significant sources of error arise 
from the calibration of the trapping power at the focus 
of the trap (7%), the asphericity of the particles (3%), 
resulting in a 6% error in the viscosity, and the uncer- 
tainty in the measurement of the particle size (1.4%), 
which contributes to a 4% error in the viscosity. The 
total error for the viscosity measurement from all error 
sources is estimated as being 10%. 
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FIG. 2: a. Signal recorded by the linear polarization mea- 
surement apparatus during rotation of a vaterite crystal. The 
frequency of rotation is (94.0 ±0.7) Hz. b. Signal traces from 
the circular polarization measurement apparatus during rota- 
tion of a vaterite crystal. 

The inferred viscosity was found to be independent of 
rotation rate, although for trapping powers above ap- 
proximately 100 mW, the viscosity was observed to de- 
crease with increasing power. This is expected, due to 
heating of the particle which results from the absorption 
of light through the sample, and also heating of the sur- 
rounding liquid due to absorption of the trapping beam 
by water, which is on the order of 10 K per W at the 
focus li^l . Figure |21 demonstrates the behaviour of the 
measured viscosity as a function of the power at the trap- 



ping focus. It was observed that higher trapping powers 
could be used with smaller particles before the onset of 
the decrease in viscosity, which is consistent with volume 
absorption effects. The observed decrease in viscosity 
with power was initially approximately linear and was 
approximately 0.13% per mW, which is equivalent to a 
temperature rise of around 0.06°C per mW. The total 
power loss for the trapping beam passing through a par- 
ticle was measured for a sample of vaterite particles, and 
was found to be within the range of 0.6% to 1.2%. The 
reflection loss due to refractive index mismatches at the 
particle surfaces is approximately 0.6%. 
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FIG. 3: Variation in viscosity with trapping power for a va- 
terite crystal of diameter 2.41 fim. The viscosity is indepen- 
dent of power for powers at the trap focus less than approxi- 
mately 100 mW. 

We have also demonstrated the ability to measure 
the viscosity inside a small confined region, such as a 
cellular membrane. Hexane filled vesicles were formed 
by emulsifying hexane, containing a small amount of 
soy lecithin (lg/1), with water containing a quantity of 
vaterite spheres. Spherical membrane structures were 
formed with diameters up to approximately 20 /xm which 
occasionally contained single vaterite crystals. The vis- 
cosity inside a 16.7 /im diameter vesicle was measured 
with a vaterite particle of diameter 3.00 /im, using the 
same procedure as outlined earlier, and was found to 
be (2.7 ± 0.5) x 10~''Pas for powers less than 100 mW, 
which is in good agreement with the established value of 
viscosity for hexane at 23°C of 3.07 x 10~^Pas (extrap- 
olated) [2^. This value is significantly less than that of 
water (9.325 x lO'^'Pas at 23°C 0. Figure H shows 
the hexane filled vesicle with the vaterite probe particle 
inside. 

The viscous drag torque acting on a spherical probe 
particle of radius a rotating with an angular frequency 
at the center of a sphere of radius R of fluid of viscosity 
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FIG. 4: Free vaterite crystal located inside a hexane-fiUed 
lipid- walled vesicle of 16.7 fim in diameter. 

/^i, surrounded by a fluid of viscosity fi2 is H^l 

TD^8nfia^nR^/{R^ -a^ + a^fii/H2)- (4) 

This can be used to determine the viscosity within a vesi- 
cle if the viscosity of the surrounding fluid is known, or 
to estimate the error if equation Q is used instead. For 
other geometries, such as if the probe particle in not in 
the centre of the vesicle, or is near a plane interface [26| . 
the effect will be more complicated. However, in most 
cases the effect of the boundary can be ignored if the 
probe particle is at least one diameter away from the in- 
terface. The effect of nearby boundaries was determined 
experimentally by measuring the viscous drag torque at 
varying distances from a solid interface. For rotation 
about an axis parallel to the interface, no change in the 
drag torque was observed until the probe particle was 
within one diameter of the interface. For rotation about 
an axis normal to the interface, no change in the drag 
torque was observed until the probe particle was very 
close to the interface. For the above measurement of the 
viscosity of hexane, the drag in the vesicle is only 0.4% 
greater than in pure hexane. The rotation rate of the 
vesicle is = + ii2{R^ - a^)/{nia^)] [H, equal 
to ri/522 for the case above. Small particles in the sur- 
rounding fluid near the vesicle were observed to orbit the 
vesicle at approximately this rate. 

The most signiflcant contribution to the error in these 
measurements is the uncertainty in accurately determin- 
ing the diameter of the spherical particle. It was not 
found to be possible to drag a particle through the mem- 
brane, and hence the method of size measurement using 
two spheres in contact could not be used. The diameter 
was instead directly estimated from an image recorded 
of the sphere. The error in determining the diameter is 
estimated as being 5%, which is the main contribution to 
the total measurement error of 19%. 



It is feasible to produce smaller vaterite particles, and 
the use of these would allow the viscosity to be probed 
on an even smaller size scale — potentially probing vol- 
umes with a capacity of a cubic micron. The ability 
to functionalise the probe particle surface would enable 
the selective attachment of the particles to various bi- 
ological structures, which would allow the torsional re- 
sponse of these structures to be investigated quantita- 
tively. We note that the torque and angular deflection 
can still be measured accurately in the absence of con- 
tinuous rotation. The demonstrated ability to measure 
viscosity within micron-sized volumes makes the system 
developed here of great value for probing the microrhe- 
ology of liquid-based materials, such as colloids. 
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